Dark soliton decay due to trap anharmonicity in atomic Bose-Einstein condensates 



o 
o 

1—5 

(N 



C/2 



I 

o 



> 
oo 

o 

ON 

o 



X 



a simple acceleration-squared law [25|, |27|, |29|, |3j 
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A number of recent experiments with nearly pure atomic Bose-Einstein condensates have con- 
firmed the predicted dark soliton oscillations when under harmonic trapping. However, a dark 
soliton propagating in an inhomogeneous condensate has also been predicted to be unstable to 
the emission of sound waves. Although harmonic trapping supports an equilibrium between the 
co-existing soliton and sound, we show that the ensuing dynamics are sensitive to trap anharmonic- 
ities. Such anharmonicities can break the soliton-sound equilibrium and lead to the net decay of the 
soliton on a considerably shorter timescale than other dissipation mechanisms. Thus, we propose 
how small realistic modifications to existing experimental set-ups could enable the experimental 
observation of this decay channel. 

PACS numbers; 03.75.Lm, 03.75.Hh, 05.45.Yv 



I. INTRODUCTION 

Dark solitons are one-dimensional phase defects that 
play an important role in nonlinear transport. They 
have been observed experimentally in nonlinear optics 
[H, shallow liquids [1], magnetic films [H and, most 
recently, ultra-cold atomic Bose-Einstein Condensates 
(BECs) H i, B 0, i, i, [M Uli- In the latter case 
matter- wave dark solitons consist of a localised dip in the 
atomic density and a phase slip, and are supported by re- 
pulsive atomic interactions. Dark solitons are prone to 
several dissipation mechanisms. Two such mechanisms 
are thermal dissipation [l^, IH, H, EBl and the snake 
instability whereby the soliton decays into vortex rings 
[S [13, 111, Ellla, Although numerous initial ex- 
periments showed such instabilities to severely limit the 
soliton lifetime [1, H, H, 0] , they can nonetheless be heav- 
ily suppressed at ultracold temperatures and in highly- 
elongated trap geometries. Indeed, such techniques have 
led to the recent generation of long-lived dark solitons 
that made many complete oscillations in the trapped con- 
densate Matter wave dark solitons are also subject 
to an instability to acceleration through an inhomoge- 
neous system [ll, [H [H [H [H . We will demonstrate 
that this instability is a generic occurrence for dark soli- 
tons in trapped atomic BECs. Furthermore, we will show 
that, under appropriately engineered conditions, this can 
become the dominant decay mechanism and is within 
reach of being detected by current experiments. 

Condensate dynamics at zero temperature obey a non- 
linear wave equation known as the Gross-Pitaevskii equa- 
tion 11]. In one-dimension and under a homogeneous 
potential, a dark soliton is an exact and integrable solu- 
tion of this equation. The presence of an inhomogeneous 
potential (which is relevant in practically all BEC experi- 
ments) simultaneously accelerates the soliton and breaks 
this integrability, rendering it unstable and leading to the 
emission of sound waves [24| . Providing the potential is 
slowly varying the decay of the soliton energy Es follows 



(1) 



where Xg is the soliton position and k is a coefficient 
that typically depends weakly on the local density and 
soliton speed. An analogous instability arises for dark 
solitons in nonlinear optics, induced by modifications in 
the optical nonlinearity, and led to the original derivation 
of the above power law [23, [1^ . 

Importantly, for the common case of harmonic trap- 
ping, dark solitons have been predicted to undergo con- 
servative dynamics in the zero temperature and one- 
dimensional limit [1^, [U [U [s^l, oscillating back and 
forth in the trap. This appears to oppose the predic- 
tion of instability under acceleration. This anomaly can 
be reconciled by the fact that in a harmonic trap con- 
tinuous reabsorption of the emitted sound stabilises the 
soliton against decay [1^. The importance of soliton- 
sound interactions was reinforced in [s^ . where it was 
shown that energy could be pumped into a dark soliton 
by driving a sound field in the condensate. It is interest- 
ing to note that a similar instability arises for a vortex 
accelerating in a trapped BEC and that they too appear 
to be stabilised within harmonic traps (ssj . 

Harmonic traps appear to be a special case. Previous 
results indicate that deviations from harmonic trapping 
tend to induce dissipation of the soliton. Such instability 
has been observed, for example, for a dark soliton in mod- 
ified harmonic traps [23l. [2511 . interacting with localised 
barriers [20l. [3^31^351. l36l. |37||. and moving through an 
optical lattice 3a 39f and disordered potential [i^l . Mo- 
tivated by this effect we study the soliton dynamics under 
several pertinent forms of confining potential. Crucially, 
we demonstrate that the presence of anharmonicities can 
destroy the soliton-sound equilibrium and lead to net de- 
cay of the soliton. This decay can be rapid and within 
current experimental soliton lifetimes, provided the gas 
is kept at a relatively low temperature. Importantly, this 
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paper demonstrates that the experimental observation of 
this decay mechanism is within reach of current long-lived 
soliton experiments [1, d, ■ 

We begin in Section II by outlining the theoretical 
framework and reviewing key properties of dark solitons 
and their dynamics in harmonic traps. Section III iden- 
tifies various experimentally-relevant geometries for ob- 
serving/inducing the emission of sound waves, firstly by 
introducing a finite cut-off to the harmonic trap, and then 
by focusing on soliton dynamics under anharmonic (gaus- 
sian; quartic) traps. In Section IV we relate our findings 
to the recent long-lived soliton experiments [sl, and 
suggest routes to detect this decay mechanism with only 
minor changes in their respective experimental set-ups. 
Our main conclusions are then summarised in Sec. V. 



numerically evaluate a meaningful local estimate as out- 
lined in Appendix A. We present our results in terms of 
harmonic oscillator units where time, l ength a nd energy 
are measured in units of ijJ~^, Iz = ^Jh/mujz and hu)z, 
respectively. The profile of the background condensate 
on which the soliton is imposed is specified by the ratio 
{fi/huJz)- For fi ^ fv-Oz, the system is in the Thomas- 
Fermi regime. Here the density profile can be approxi- 
mated by n(r) — — V{r)]/g and the system is much 
wider than the soliton. In the opposite limit of /i <C hujz 
the system is closer to the non-interacting limit with size 
comparable to the soliton size. Throughout this work and 
unless otherwise stated we will employ a system defined 
hy fi — 22. AhiUz, which is well within the Thomas-Fermi 
regime. 



II. BACKGROUND & MODELING 



B. Properties of dark solitons 



A. Gross-Pitaevskii equation 

At ultra-cold temperature and under weak atomic s- 
wave interactions the vast majority of the atoms in an 
atomic BEG reside in the same quantum state. The 
atomic BEG can then be parameterised by a macro- 
scopic 'wavefunction' ip{r,t). This can be expressed as 
■^(r, t) = ■\/7i(r, t) exp[i6{r, t)], where n(r, t) is the atomic 
density distribution and 9{r,t) is the condensate phase 
which can, in turn, be related to a fluid velocity via 
v{r,t) = {h/m)Wd{r,t). This wavefunction is described 
by the Gross-Pitaevskii equation (GPE) [4l| . 
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dt 



= [ -^V2 + K,t(r) + ,g|V(r, t)p ) i,{v, t).{2) 



Here m is the atomic mass and g — 47rfi^ ag/m charac- 
terises the nonlinear atomic interactions, where is the 
s-wave scattering length. The chemical potential of the 
condensate is denoted by /x. The external confining po- 
tential Vcxt(r) is typically harmonic with the form, 
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In ID and on a uniform background density rtg a dark 
soliton with speed v and position (z — vt) has the ana- 
lytical form. 



^{z,t) = v^e-'(^/'''* ^/3tanh 



/3 



(z - vt) 



+ i 



where /3 = -^/l — (w/c)^ and the healing length ^ = 
h/ y/mnog characterises the soliton width. For this homo- 
geneous case the chemical potential is given by = nog. 
The soliton speed v is related to both the soliton depth 
rill (with respect to the background density) and the to- 
tal phase slip S across the centre via v — \/no — rid = 
ccos(S'/2), where c = \friQgJrn is the Bogoliubov speed 
of sound. A stationary soliton has a tt phase slip and a 
maximum depth = ng, while a soliton with speed c is 
indistinguishable from the background fluid. 

The energy of the dark soliton of Eq. (jl]) , which is an 
ideal quantity for parameterising the soliton behavior, is 
given by. 
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where we have assumed cylindrical symmetry, and where 
LOr and ujz are the harmonic trap frequencies in the radial 
and axial directions, respectively. Assuming a highly- 
elongated condensate (uJz <C cj^) with tight radial con- 
finement (hiOr > m)i the condensate dynamics become 
quasi-one-dimensional and the ID version of the GPE 
[H, [H, H, [H becomes applicable. (This also assumes 
that fluctuations in the condensate phase can be ignored, 
which is justified at such low temperatures [13].) 

We will simulate the condensate dynamics using the 
ID GPE. This is performed by numerically propagating 
the wavefunction ^{z,t) on a discrete spatial grid using 
the Grank-Nicholson technique [46]. Although the en- 
ergy of the soliton technically extends to infinity, we can 



Accordingly, the dissipation of energy from the soliton is 
associated with an increase in its speed. Under an inho- 
mogeneous potential the soliton behaves, to first order, as 
an effective particle of negative mass [1^ . The soliton will 
tend to oscillate when confined within a trap and any dis- 
sipation is associated with an increase in the amplitude of 
the soliton oscillations, an effect termed "anti-damping" 
[2^ . For the specific case of a harmonic trap a dark soli- 
ton ispredicted to oscillate at a frequency ujs = i^z/V2 
[H, m m, Hi m, m, [M HI . Note that deviations from 
this prediction can arise due to the dimensionality of the 
system j43| , the presence of other solitons and trap anhar- 
monicities, and these effects have been recently observed 
experimentally [1, . 
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FIG. 1: (Color online) Dark soliton dynamics in a 
harmonically-trapped condensate for a soliton initially lo- 
cated at the trap centre and with a speed vo = 0.5c. (a) 
Renormalised density |?/)(2;,t)p — nTi(z), in units of the peak 
density no. (b) Soliton position Zs{t) (solid line, left axis) and 
condensate centre of mass Zcom(i) (dashed line, right axis). 
The soliton position Zs is defined as the point of the soliton 
density minimum, (c) Soliton energy Ea{t), determined nu- 
merically as outlined in Appendix A. 



C. Dark solitons in purely harmonic traps 



We first review the dynamics of a dark soliton in a ID 
harmonic trap V{z) = ma;^z^/2, which has been the sub- 
ject of much theoretical work 0, 

The initial state is ip{z,0) = ^/nTI(z)'4>siz,0), where 
nTi(z) is the time-independent background density ob- 
tained by imaginary time propagation of the GPE [44] 
and the soliton solution ipsiz, 0) is given by Eq. (|4]). The 
dark soliton solution is centred at the origin, where the 
density is locally homogeneous, and has initial speed vq. 

The typical evolution of a dark soliton in a 
harmonically-confined BEG is shown in Fig. [TJa). The 
soliton (dark density dip) oscillates sinusoidally through 
the BEG. This is accompanied by small-scale density 
excitations of the background condensate, of amplitude 
^ O.lno. We attribute these to the sound waves emitted 
from the accelerating soliton. The sound waves remain 
confined within the system and form a dipole (sloshing) 
mode of the condensate. We have tracked the soliton tra- 
jectory Zg (solid line in Fig. [TJb)) and taken its Fourier 
spectrum (solid line in Fig. [^J^a)). This confirms that the 
soliton primarily oscillates at frequency ^ '^/■/^/^^ ^-^ 
agreement with previous predictions [12, [23, |3i| . 
The spectrum features a secondary frequency component 
at Lu = ujz, which arises due to the interaction of the soli- 
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FIG. 2; (a) Fourier spectrum of Zs{t) (solid line) and Es{t) 
(dotted line) for the case in Fig. [1] (b) Average soliton fre- 
quency ujs as a function of the system parameter {fujz/l^) for 
Vq — 0.25c (squares), vq — 0.5c (crosses) and vq = 0.75c (cir- 
cles). The analytic prediction, uj^^ = ujz/\/2, is shown (dashed 
line) . 



ton with the dipole mode of the condensate. 

During these dynamics the condensate centre-of-mass 
oscillates due to the linear momentum introduced by the 
imposed soliton solution. The condensate centre-of-mass 

^com 

= z\^{z)\Mz/N is plotted in Fig.^b) (dashed 
line, right axis). The centre-of-mass motion is decoupled 
from the "internal" dynamics of the condensate and os- 
cillates sinusoidally at the trap frequency ujz, as noted in 
(23I I. The soliton energy Eg (sohd line in Fig.[Hc)), whose 
definition is outlined in Appendix A, remains constant 
on average, but features oscillations. The corresponding 
Fourier spectrum (dotted hue in Fig.[2ja)) has contribu- 
tions from 2lUs and beat frequencies between the soliton 
and dipole mode, e.g. {lOz — ajg), {iVz + i^s), plus higher 
order combinations (not visible in Fig. [2]). 

The system considered here is in the Thomas-Fermi 
regime /i ^ huz (with /i = 22Ahujz)- We have also anal- 
ysed the soliton motion over a range of values of {fx/hiUz) 
with the results shown in Fig. [DJb). For /i 3> hujz the 
soliton oscillation frequency lOs agrees with the analytic 
prediction of Ws ~ iOz /V2 12, 16,,^,^]. Indeed, these 
predictions are based on treating the soliton as a clas- 
sical particle and so is valid when the soliton is small 
compared to the system size. As the system parameter 
(fi/hujz) is reduced (which reduces the size of the system 
relative to the soliton size), lOs deviates from this predic- 
tion and becomes dependent on the soliton speed. For 
fi = 5hujz, which is when the system size is comparable 
to the soliton size, uj^ can be over 10% larger than the 
analytic predictions. Note that in all cases considered, 
the soliton dynamics undergo no net decay. Furthermore, 
the frequency component at uJz is present throughout the 
range of systems tested. Note that shifts in the soliton 
oscillation frequency can also arise due to the presence of 
the transverse direction [1, [4^ and this is distinct to our 
present ID analysis. 

We attribute the density excitations in the system to 
sound waves emitted by the unstable accelerating soliton. 
Although these waves co-exist with the soliton at long 
times, it takes a finite time from t — for the emitted 



sound waves to reflect ofl' the edge of the condensate and 
reinteract with the soHton. This will be approximately 
half a trap period t « ttwJ^. Close inspection of Fig.[Uc) 
shows that the soliton energy decays monotonically at 
these early times. Indeed, for this time the decay follows 
the acceleration-squared law of Eq. ([T]). After this time 
the soliton energy is found to oscillate between its initial 
value and some lower value, providing a clear indication 
of the reabsorption of sound by the soliton. 

A recent study by Pelinovsky et al. [1^1 has consid- 
ered the motion of a dark matter-wave soliton in a har- 
monic trap by means of an asymptotic multiscale expan- 
sion method. Sound radiation in the system is predicted 
to oscillate at the same frequency as the soliton, thereby 
causing the periodic re-interaction of the sound and soli- 
ton, and thus enabling reabsorption by the soliton. This 
prediction is somewhat consistent with our observations, 
although we observe the sound waves to have oscillatory 
components at the trap frequency as well as the soliton 
frequency. 

We now discuss various geometries in which the emit- 
ted sound waves can be studied. 



III. PROBING ANHARMONICITY-INDUCED 
SOLITON DECAY 

A. Dark soliton dynamics in a cut-off harmonic 
trap 

Although a purely harmonic trap generates no net de- 
cay of the dark soliton, we can nonetheless probe the 
existence of sound emission even within harmonic traps 
by introducing a finite cut-off to the trap, as discussed 
in [2^ . We thus consider a cut-off harmonic trap of the 
form. 



(a) 10 



for 



\z\ < y /2Vo/mujl 
\z\ > ^2Vo/mcjl 



(6) 



where Vq specifies the magnitude of the cut-off. Although 
this potential is somewhat unphysical, it can be approxi- 
mated by employing a tight dimple trap embedded within 
a weaker ambient trap (i.e. an effective double-trap geom- 
etry) [25I I . Since sound waves have an energy of approxi- 
mately the chemical potential /i, the cut-off depth enables 
control between a regime where the sound waves remain 
confined to the inner trap (Vq 3> fJ-) and a regime where 
they can escape the inner trap (Vq < fJ.). In Fig. [3] we 
show the evolution of the system for various values of Vq 
which illustrate these regimes. For Vq = 1.2^ (Fig. [3]^ a)) 
the sound waves are confined to the harmonic region and 
the evolution is essentially identical to that under pure 
harmonic trapping (Fig. [TJa)). For this case the soliton 
energy, presented in Fig. Ul^a), shows no net decay. For 
Vq = 0.8/x (Fig. [3{c)) all of the emitted sound waves es- 
cape the harmonic region and propagate outwards. The 
soliton energy (Fig. IUJa)) decays monotonically and fol- 
lows the acceleration-squared law of Eq. ^ [2^. Note 
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FIG. 3: (Color online) Space-time plots of renormalised den- 
sity {\Tp{z,t)\^ — nTi{z)) showing the dynamics of a = 0.5c 
dark soliton in cut-off harmonic trap with (a) Vo/fJ^~ 1.2, (b) 
1.05 and (c) 0.8. The dark soliton appears as a black line, 
while sound waves appear as light/dark areas. 



that at late times the soliton has anti-damped sufficiently 
that it escapes to the outer homogeneous region where 
it propagates with constant speed. In between these 
regimes, e.g. for Vq = 1.05^ (Fig. [3Kb)), only a fraction 
of the sound waves escape to infinity, and we observe a 
slower decay of the soliton. Note that the power law of 
Eq. ([1]) no longer applies since considerable sound re- 
absorption occurs. The emitted sound remains mostly 
confined to the central trap and the amplification of this 
sound field as the soliton decays is clearly evident. 

For the cut-off harmonic trap, decay of the soliton oc- 
curs because sound waves escape the harmonic region 
and carry away their energy. For Vq < this can al- 
ways happen. However, it is important to understand 
how decay changes as Vq becomes greater than as this 
will prove crucial in our subsequent discussion on soli- 
ton decay in experimentally relevant gaussian traps. We 
define a decay time td as the time for the soliton to be- 
come indistinguishable from the background condensate. 
In Fig. Hfb) we show how the decay rate f^^ varies as 
Vq becomes greater than fj,. The most important thing 
to note is that the decay rate drops very quickly with 
Vq. For the soliton speed vq = 0.5c, which we have con- 
sidered in Fig. [31 decay ceases for Vq > 1.15/i. This 
corresponds to the point at which sound waves no longer 
escape the harmonic region. This threshold is larger for 
faster solitons. For example, for vq = 0.75c, decay is 
ultimately prevented for Vq > lAfi. This suggests that 
the sound waves emitted by faster solitons are of higher 
energy. However there always exists a threshold cut-off 
depth beyond which the tendency to decay is suppressed. 

Having reviewed and extended previous results of dark 
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FIG. 4: (a) Evolution of the soliton energy Es in a cut-off 
harmonic trap for Vo/fi = 0.8, 1.05 and 1.2 for initial speed 
vo = 0.5c. (b) Decay rate as a function of trap cut-off Vo 
for Vo = 0.25c (circles), vo = 0.5c (crosses) and vo = 0.75c 
(pluses) . 



soliton dynamics in harmonic traps we now progress to 
consider an important additional decay mechanism for 
dark solitons, namely their propagation in condensates 
under anharmonic confinement. 



B. Dark soliton dynamics in a gaussian trap 
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FIG. 5: Space-time plots of renormalised density {\'ip{z,t)\^ — 
nTi{z)) showing the dynamics of a t;o = 0.5c dark soliton in 
a (a) deep gaussian trap with Vo = 5/i, (b) shallow gaussian 
trap with Vo = 2/i, and (c) infinite system with Vo = 0.8/i. 



For simplicity, and in order to make contact with re- 
alistic experiments, we begin by considering a gaussian 
trap; such traps are commonly used experimentally and 
are formed in the cross-section of an off-resonant laser 
beam. Notabl}^j_^the recent dark soliton experiment of 
Becker et al. &] was performed in a trap which was 
gaussian in the longitudinal direction. We will focus on 
this particular experiment further in Section V. However, 
in this section we will present the general dynamics in a 
gaussian trap. 

We consider here a ID gaussian trap of the form. 
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The leading term in the Taylor expansion is harmonic, 
with trap frequency lOz, while the higher order terms are 
anharmonic. In a similar manner to the cut-off trap, 
the parameter Vq controls whether sound waves can es- 
cape the system. However, we will generally consider the 
regime Vo > ^J.^ where the condensate (including soliton 
and sound) is wholly confined to the gaussian trap. Im- 
portantly, the gaussian depth Vq also controls the degree 
of anharmonicity experienced by the condensate, and can 
be varied experimentally by varying the laser intensity. 
In the limit Vb 3> /i, the confinement is effectively har- 
monic, while for Vb ~ /i the condensate will experience 
strong anharmonic terms. 

We first consider how the gaussian depth affects the 
soliton dynamics, with Fig. [5] showing the density evolu- 
tion for three values of Vq. For Vb = 5/i [Fig. (Ha)], the 



trapping is effectively harmonic and the dynamics are 
similar to those in a harmonic trap, with the soliton en- 
ergy, shown in Fig.[B]Ja), undergoing no net decay. As for 
the purely harmonic trap this occurs due to the complete 
reabsorption of the emitted sound by the soliton. 

In contrast, for Vq = 2/i (Fig. (H^b)), the soliton am- 
plitude increases and its energy decreases (Fig. El^a)). 
The soliton decay is accompanied by the growth of sound 
waves in the system. After only 5 oscillations {t ~ 50wJ^) 
the soliton is so fast and shallow that it is indistinguish- 
able from the sound field. It is important to note here 
that the trap depth is sufficiently large that no sound 
waves escape the system and so the decay mechanism is 
quite distinct from that of the shallow cut-off harmonic 
trap. While the sound waves remain confined to the trap 
and co-exist with the soliton, they do not stabilise the 
decay. We believe that the anharmonic nature of the 
trap introduces a dephasing between the sound and soli- 
ton that prevents the sound from being completely reab- 
sorbed. For completeness we also present the dynamics 
for Vo = O.Sfi (Fig. [DJc)) where sound waves can escape 
the trap. In this case we observe the soliton undergoing 
decay with no sound reabsorption, similar to the corre- 
sponding case for the cut-off harmonic trap (Fig. [3jc)). 

In Fig. ^h) we plot the decay rate ij^ of the soliton 
versus the gaussian depth Vq for various soliton speeds. 
Consider the case of wo — 0.5c (crosses). Here decay 
occurs up to Vo = 3.5fi. Note that according to the 
results of Fig.lDJb) we expect that no sound can escape for 
Vo > l-5fi and so it is clear that the decay occurs even in 
the presence of the emitted sound. We see that the decay 
rate decreases with Vb . This is consistent with the decay 
being due to anharmonicities: as Vb increases the trap 
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FIG. 6; (a) Evolution of the soliton energy Es in a gaussian 
trap for various depths Vb for a soliton with initial speed vq — 
0.5c. (b) Decay rate as a function of Vb for vo = 0.25c 
(circles), vq — 0.5c (crosses) and «o ~ 0.75c (pluses). 
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FIG. 7; (a) Evolution of the soliton energy Es in a quartic 
trap for soliton speeds of (from top to bottom) vq/c — 0.5, 0.6, 
0.7 and 0.8. The trap has the form of Eq. ^ with Vb = 10^. 



becomes increasingly harmonic (from the perspective of 
the condensate). Furthermore we see that faster solitons 
have larger decay rates and undergo decay up to larger 
values of Vq. We attribute this to the fact that faster 
solitons probe closer to the edge of the condensate where 
the trap anharmonicities have the largest effect. 

It is important to note that the soliton decay in the 
gaussian trap can be considerable. For example, for the 
case presented in Fig. O^b), the soliton has practically 
disappeared hy t = SOwJ^, which corresponds to approx- 
imately 5 soliton oscillations in the trap. Assuming a typ- 
ical trap frequency of oj^ = 27r x 100 Hz, this timescale 
corresponds to around 80 ms, which is well within the 
thermodynamic lifetime of dark solitons at ultralow tem- 
peratures fs", TB*! . Furthermore, the fastest decay pre- 
sented in Fig. [BJb) occurs twice as fast as this example, 
corresponding to only two or three soliton oscillations. 

C. Dark Soliton Dynamics in a Quartic Trap 

To give further evidence of the role of trap anharmonic- 
ities in dark soliton decay we briefly consider the soliton 
dynamics within a quartic trap. Quartic potentials can 
be obtained, for instance, by modifying the gaussian po- 
tential from a laser beam, as performed experimentally 
to create harmonic-plus-quartic traps [IJ. As such, we 
assume our quartic trap to have the form of the quartic 
term of the gaussian potential (Eq. [5]), i.e.. 

Following previous sections, we define our system via 
fi = 22AfujJz, and assume here Vq = 10/x. In Fig. [7] 
we present the evolution of the soliton energy for vari- 
ous soliton speeds. We observe that the soliton stability 
is strongly dependent on speed. Slow solitons vq < 0.5c 
undergo little or no detectable decay, while faster soli- 
tons, e.g. Vo — 0.8c, do undergo decay. Note that the 
decay is significantly slower than in the gaussian trap. 

We have also considered further traps formed by a 
combination of harmonic and quartic terms, and observe 



regimes of soliton decay in all cases. For a harmonic- 
plus-quartic trap, the trap has similar properties to the 
purely quartic trap and we obtain similar qualitative re- 
sults. For a harmonic-minus-quartic trap, the trap is 
similar to a gaussian trap, particularly in the sense that 
it only has finite depth, and we find similar results to the 
gaussian trap. 

We now discuss how our findings relate to recent ex- 
periments conducted at very low temperatures. 

IV. APPLICATION TO RECENT 
EXPERIMENTS 

In the early dark soliton experiments 0, @, the 
solitons were short lived due to thermal dissipation and 
the snake instability. This limited the observation times 
to a fraction of the soliton oscillation period. Recent 
experiments [1, [§, have heavily suppressed these dis- 
sipation mechanisms by operating at extremely cold tem- 
peratures and by employing quasi-lD geometries. As a 
result they were able to observe solitons over much longer 
timescales in which the soliton undergoes multiple oscil- 
lations in the trap. Indeed, in the experiment of Becker 
et al. ^ dark solitons were observed for over 2.8 s, corre- 
sponding to approximately 10 soliton oscillations in the 
trap. The considerable soliton lifetimes in these recent 
experiments, as well as the suppression of thermal dissi- 
pation and the snake instability, offer realistic opportu- 
nities to experimentally detect the acceleration-induced 
decay of the soliton discussed here. 

A. Experiment of Becker et al. 

We firstly consider the experiment of Becker et al. [§] . 
Dark solitons were generated in an elongated ^^Rb con- 
densate via a phase imprinting technique and observed to 
oscillate in the trap with peak velocity wq — 0.56c. The 
trap, formed by the optical dipole force of an applied 
laser beam, featured a gaussian profile in the axial direc- 
tion with waist Wz = 125 /xm and dominant harmonic 
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frequency = 27r x 5.9 Hz. For the experimentally- 
quoted peak density of uq = 5.8 x 10^^ cm"^ we es- 
timate the chemical potential of the condensate to be 
fi/ks ^ n^g/ks ^ 23 nk. The amplitude of the gaus- 
sian potential is Vb = muilW^/A — 58 nK. The impor- 
tant depth-to-chemical potential ratio is then Vq/ fi ^ 2.5. 
Based on these parameters our simulation gives a soliton 
oscillation frequency of = 27r x 3.94 Hz and an ampli- 
tude of Zq « 34/xm, which are both in good agreement 
with the values quoted in Ref . . 

It is pertinent to ask whether acceleration-induced de- 
cay of the soliton occurs in this system, and whether it 
may be detectable. In Fig. EJa) we plot the soliton tra- 
jectory in the experimental system for several values of 
the trap depth Vq. Note that the trap depth Vq is propor- 
tional to, and can be varied by modifying, the intensity 
of the laser beam that forms the trap. For Vo = 4/i 
(grey line) the trap is essentially harmonic: the soliton 
undergoes no discernable decay and oscillates with con- 
stant amplitude. For Vq = 2^ (dotted hne), which is just 
smaller than the number used in the experiments, the 
soliton experiences weak trap anharmonicities, causing it 
to decay slowly (as evident from the increase in the soli- 
ton amplitude) and effectively disappear by i ~ 6 s. This 
decay timescale is longer than the reported experimental 
lifetime (2.8 s), and so evidence of such decay would not 
be easily visible in the current experiment. Under a slight 
reduction in the trap depth to Vq/^ = 1.7 (solid line) the 
soliton undergoes much quicker decay, disappearing into 
the background condensate after only 2.5 s. Importantly, 
this is now within the experimental soliton lifetime. 

The effect of trap depth is shown in more detail in 
Fig. ^h) which plots the decay time as a function of 
trap depth Vo for two different initial soliton speeds. The 
maximum reported experimental soliton lifetime (2.8 s) 
is shown (dashed line) for reference. As earlier results in- 
dicated, the decay time generally increases with Vq as the 
trap becomes increasingly harmonic. For the experimen- 
tal trap depth Vq ~ 2.5fi and soliton speed vq = 0.56c 
the soliton decays over a timescale id ~ 6 s, which is 
longer than the experimental lifetime. How could one 
modify the experiment to make the decay occur within 
the experimental timescale? Firstly, this can be achieved 
(assuming an initial soliton speed of vq = 0.56c) by re- 
ducing the trap depth to the range Vq ^ 1.8^. Secondly, 
the decay can be made more rapid by beginning with a 
faster soliton. For example, for vq = 0.75c the soliton 
will decay within the experimental lifetime provided the 
trap depth does not exceed Vq ^ 2.8/x. 

These results suggest that with only slight modifica- 
tions this decay channel will dominate the experiment 
and cause the soliton to decay within the experimental 
timescale. The soliton dynamics are also very sensitive 
to the trap depth Vq, a feature which also lends itself 
to experimental detection. The detectable signature for 
the decay would be provided by the considerable anti- 
damping of the soliton that accompanies the decay, which 
could be probed either in expansion images [1, d, [^, 
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FIG. 8: Soliton decay in the experimental system of Ref. @. 
(a) The evolution of the soliton position for Vo = 1.7fi (solid 
line), Vo = 2/i (dotted line) and Vo — ifi (grey line). Following 
the experimental observations we employ an initial soliton 
speed of vo = 0.56c. (b) Decay time in the experimental 
system as a function of the gaussian trap depth Vq. Initial 
soliton speeds of vq = 0.56c (circles) and vq — 0.75c (crosses) 
are presented. The dashed line indicates the maximum soliton 
lifetime reported in the experiment. 



or even in situ 



B. Experiment of Weller et al. 

We will also briefly discuss the recent soliton experi- 
ment by Weller et al. Q. Two separated condensates 
were initially formed in a double trap potential. The 
double trap potential is formed by the combination of 
an elongated harmonic trap and one period of a one- 
dimensional optical lattice potential. Upon removal of 
the optical lattice the condensate underwent nonlinear 
interference which resulted in the production of dark soli- 
tons. These solitons were then observed to make up to 7 
oscillations in the system. Although the solitons are typ- 
ically produced in pairs we expect our predictions to hold 
for multiple solitons. The final trap potential is harmonic 
and so we do not expect any net decay of the solitons due 
to acceleration-induced sound emission. However, if the 
optical lattice potential were instead maintained at some 
low amplitude following the creation of the dark solitons, 
then we anticipate that soliton decay would occur. 

We previously considered the dynamics of a dark soli- 
ton in a harmonic trap featuring an optical lattice po- 
tential [111. The optical lattice potential was assumed 
to be weak (amplitude less than the chemical potential) 
such that the soliton could traverse the lattice site/s 
and thereby undergo acceleration-induced sound emis- 
sion. Provided the lattice spacing A was in the range 
^ < A < 2Lz, where ^ is the healing length and is 
the condensate length, the soliton consistently underwent 
decay. (For A < ^ the density perturbations due to the 
optical lattice become smeared out while for A > 2Lz 
the trap reduces to being effectively harmonic.) This is 
consistent with the results presented in this paper: the 
optical lattice introduces anharmonicities that prevent 
complete sound reabsorption. The ensuing decay could 
be rapid, with complete decay occurring as quickly as 
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within one trap period in certain cases. In the experiment 
of WeUer et al. [1] the lattice spacing is approximately 
Lz , and so we expect that the presence of the lattice can 
indeed induce soliton decay. This suggestion is similar 
to an earlier prediction (20| where we showed that a sin- 
gle gaussian barrier within a harmonic trap could lead to 
soliton decay. Moreover, decay in the experiment of [sj 
could be further enhanced by decreasing the periodicity 
of the proposed optical lattice; following previous results 
[11] we would expect the quickest decay when A ~ 5^. 



V. DISCUSSION 

We have examined the zero-temperature, mean-field 
dynamics of dark solitons in atomic Bose-Einstein con- 
densates under harmonic and anharmonic trapping. Such 
inhomogeneous potentials generally accelerate the soliton 
and induce an instability whereby the soliton radiates 
sound waves. In the absence of sound reabsorption the 
power radiated by the soliton is proportional to its accel- 
eration squared. However, due to the finite-size of BECs, 
the sound waves remain in the system and generally re- 
interact with the soliton. The ensuing dynamics between 
the soliton and sound waves are crucially dependent on 
the form of the trapping potential. 

Under harmonic trapping no net decay of the soliton is 
observed, despite the local instability of the soliton. The 
isochronicity of the harmonic trap leads to the creation 
of a dynamical equilibrium between the soliton and the 
emitted sound waves which prevents net soliton decay. 
We note that the inclusion of an additional harmonic fre- 
quency to the system, e.g., by considering an asymmetric 
harmonic trap (different trap frequencies for z > and 
z < 0), maintains the stability of the soliton. Decay of 
the soliton in the harmonic trap can be induced, but in- 
volves removal of some of the emitted sound, e.g. by 
introducing a finite cut-off to the harmonic trap. 

In a gaussian trap decay can occur without the removal 
of the emitted sound even at zero temperature. This is 
possible because the soliton-sound equilibrium is broken 
by the anharmonicity of the trap. In accord with this the 
decay is strongly dependent on the anharmonicity expe- 
rienced by the soliton. The soliton decay is accompanied 
by the growth of short wavelengths in the sound field. 
Quartic trapping also leads to soliton decay. Further- 
more, several existing studies have added some form of 
anharmonic perturbation to the harmonic trap and ob- 
served soliton decay, e.g. under the addition of an optical 
lattice [H, [3^, a localised barrier [20|, [3l| and a disor- 
dered potential (40| . in concurrence with our findings. 

Soliton decay can be rapid. For example, in the gaus- 
sian trap we have observed a soliton (with initial speed 
vq = 0.5c, where c is the speed of sound) disappear into 
the sound field after only five oscillations in the trap. 
Given that solitons have recently been stabilised for over 
7 trap oscillations [1, Q it is clear that acceleration- 
induced soliton decay can dominate over the more well- 



known dissipation mechanisms (thermal dissipation and 
the snake instability). For the experiment of Becker et 
al. ■9'], which was conducted in a gaussian trap, our re- 
sults suggest that the soliton should completely decay 
within the experimental lifetime provided a slightly shal- 
lower gaussian trap or a faster soliton is employed. Note 
that the trap depth is readily varied through the inten- 
sity of the laser beam that forms the trap, while the soli- 
ton speed can be varied by adjusting the degree of phase 
imprinting. Furthermore, for the recent experiment of 
WcUcr et al. [1| we suggest that if a weak optical lattice 
is maintained during soliton propagation it may lead to 
significant soliton decay. Soliton decay can be detected 
by the increase in the soliton oscillation amplitude mea- 
sured either in successive absorption images or in situ, 
with the accompanying growth of the sound field provid- 
ing further indication of this effect. 

It is interesting to note that the instability to accel- 
eration that underpins this decay channel has also been 
predicted for vortical excitations in atomic BECs, includ- 
ing single vortices [H,!!!,!!^, corotating vortex pairs [s^l 
and two-dimensional vortex-antivortex pairs |511] . Hybrid 
dark soliton/ vortex rings excitations that have recently 
been observed 152] may also be expected to undergo 
an instability to acceleration. As such the decay channel 
discussed here may extend well beyond dark solitons. 
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APPENDIX A: EVALUATION OF THE DARK 
SOLITON ENERGY 

Although the soliton energy technically extends up to 
the boundary of the system it is useful to have a local 
measurement of the soliton energy. We perform this by 
numerically evaluating the energy in a local region about 
the soliton centre. In trapped systems this will typically 
contain contributions from other excitations, e.g. sound 
waves. However, this contribution is typically very small. 
Note that from an experimental point of view it is pre- 
cisely the behaviour of this soliton region, including local 
excitations, which is probed experimentally. 

The energy density of the condensate is given by the 
energy functional, 

em^^m\' + V,^,\M' + ^\^f. (Al) 
We define our soliton energy to be, 

Es^ e[i;]dz~ eJV^dz. (A2) 

J Zg—Zij^t J ^g — 2int 
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Here z-mt defines the extent of the integration region. 
Note that we subtract ofi' the corresponding energy of 
the condensate background density (second term). Fig- 
ures [HIa) and (b) present the integrated sohton energy 
as a function of Zjnt for various sohton speeds and back- 
ground densities, respectively. For simphcity we consider 
an infinite system with V{z) = 0. For large enough Zi^t 
we observe convergence to the analytic soliton energy of 
Eq. [5] (dashed lines). Larger Zint is required for conver- 
gence for faster solitons and/or lower densities, due to 
the increased soliton width. However, for the speeds and 
densities considered here, the integrated energy is vir- 
tually indistinguishable from the analytic prediction for 
^.int _ which is employed in this work. 
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FIG. 9: Integrated local soliton energy Es as a function of 
integration width z'"'. (a) Various soliton speeds are consid- 
ered in a system with peak density n = no- (b) Various values 
of background density n are considered for fixed soliton speed 
V = 0.5c. In each case the analytic soliton energy is indicated 
by a dashed line. This is performed in an infinite V{z) = 
system and our units of energy and distance are the chemical 
potential jj, = nog and healing length ^ = h^mnog of the 
n = no system. 
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